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TECHNICAL NOTE 2132 


THE CALCULATION OF MODES AND FREQUENCIES OF A MODIFIED 
STRUCTURE FROM THOSE OF THE UNMODIFIED STRUCTURE 
By Edwin T. Kruszewski and John C. Houbolt 


SUMMARY 


A method is developed for the calculation of the natural coupled or 
uncoupled frequencies and modes of a structure with modifications, such 
as the addition of concentrated masses or springs, directly from the 
known modes and frequencies of the unmodified structure. The modes of 
the modified structure are expanded in terms of the modes of the unmodi- 
fied structure. A characteristic equation and a frequency determinant, 
the order of which is twice the number of modifications, are derived by 
the use of the Galerkin method. Numerical examples are presented to show 
the accuracy of the method and the number of modes and frequencies of the 
unmodified structure necessary for agreement with exact solutions. 


INTROnJCTION 


The calculation of the natural modes and frequencies of an airplane 
structure is usually required for various loading conditions. The varia- 
tion of conditions may be brought about by changes in pay load, changes 
in the amount of fuel carried, the addition of tip tanks, and so forth. 
These changes may be regarded as modifications to a primary structure. 

In addition to weight changes, the addition of elastic restraints, such 
as spring supports which may be used in ground vibration tests, may also 
be considered as modifications. Thus, in this paper the basic or primary 
structure is known as unmodified, and the structure after masses and 
springs are added is known as modified. 

Present methods of calculating modes and frequencies require a sepa- 
rate and independent calculation for each modification to the primary 
structure. In order to simplify the calculations of these modes and fre- 
quencies, a method is developed in this paper that allows, with very 
little extra work, the calculation of the modes and frequencies of a modi- 
fied structure directly from the modes and frequencies of the unmodified 
structure . 
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The paper presents "both a theoretical analysis and numerical 
examples . In the theoretical analysis a frequency determinant and the 
modal functions are derived. The order' of the frequency determinant is, 
in general, twice the number of modifications to the primary structure. 
Each element of the determinant is a series having terms that are func- 
tions of the modal shapes and frequencies of the primary structure. The 
numerical examples illustrate the convergence of the series and the accu- 
racy with which the frequencies and modes of the modified structure can 
he calculated. 


SYMBOLS 


El 

GJ 

y 

<p 




<Pn 

(JQ 

% 

Mi 

a i 


e i 


k 


flexural stiffness 
torsional stiffness 

deflection of elastic axis of modified "beam 

angle of' twist of cross section of modified "beam 

deflection of elastic axis of unmodified "beam in nth coupled 
mode normalized to give unit tip deflection 

angle of twist of cross section of unmodified "beam in 
nth coupled mode 

a natural frequency of modified beam 

natural frequency of unmodified beam in nth coupled mode 

magnitude of concentrated mass at station x-^ 

spring constant of spring at station x^ 

proportion of nth. mode of unmodified beam present in a 
modified-beam mode 

distance between center of gravity of cross section of beam 
and elastic axis; positive when center of gravity lies 
forward of elastic axis 

distance between elastic axis and center of gravity of con- 
centrated mass or distance between elastic axis and spring 
(see fig. 2); positive when center of gravity of mass or 
spring lies forward of elastic axis 

radius of gyration of cross section of beam about elastic axis 
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k i 


m 


Hi 

BGt 



for a mass, radius of gyration of mass about elastic axis; 
a spring, e-^ 

mass of unmodified beam per unit length 
spanwise coordinate measured from center line 
lateral load caused by mass or spring at station x^ 

torque about elastic axis caused by mass or spring at 
station x^ 

= K i ( e i Mi + k i 2 &n]i) 
f m(y n 2 + 2ey n cp n + k^/^dx 


for 


L 


semispan of beam 


Ki 


p 

= M-^co for concentrated mass at station x^ 
= -a^ for spring at station x^ 


f frequency coefficient 


8(x - XjJ function such that its value is zero at every point except 
argument x^ and its value at this point is infinity in 


such a way that 



x^)dx = 1 


Subscripts : 

i, j integers referring to locations of concentrated masses or springs 

m, n integers referring to modes of unmodified beam 

A, B, C locations of masses in example 2 


THEORETICAL ANALYSIS 


Development of characteristic equation .- Figure 1 illustrates a 
typical structure that is to be considered in the development of the 
method presented in this paper. It consists of a primary beam of any 
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spanwise -variation of flexural stiffness EE, torsional stiffness GJ, 
and mass m. The modifications to the primary structure are the addi- 
tion of any number of concentrated masses of magnitude M-^, Mt>, . . . 

at stations x^_, X 2 , . . . , respectively, and any number of springs with 
spring constants a^, ag, ... at stations x^, xg, . . . , respectively. 

The method developed in this paper is based on the assumption that 
the coupled modes y n and cp n , normalized to give a unit tip bending 

deflection, are known for the primary structure. If the uncoupled modes 
of the modified structure are desired, uncoupled modes of the primary 
structure must of course be used. 

The procedure presented in this paper is based upon the Galerkin 
method (reference l). The differential equations for the modified beam 
in free harmonic coupled vibration, vibrating at a natural, frequency cn, 
are 


[P] ± &(x - Xj) = 0 ' (la) 


[T]i&(x - Xi) = 0 (lb) 


dx 


j^ 1 - ma£(y + e<p) - ^ 


and 


g|) + ™^(ey + k3p) + 


where y and cp are the maximum coupled bending and torsional deflec- 
tions of the modified beam, respectively, and [P] i and C-^Ili are the 

concentrated loads and torques, respectively, caused by the modifications 
at station x-^ and are defined as 


Mi - K i(Edi + e i&]±) 
Hi - K i( e iEii + k 2 &]i) 


where 

r 




-a ± 


for concentrated masses 
for springs 


(i = 1, 2, 



( 2 ) 
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e.^ distance from the elastic axis to the center of gravity of 

concentrated mass or, in the case of a spring, to the center 
line of the spring (see fig. 2) 


ki 


radius of gyration of mass about elastic axis for concentrated 
masses, or, in the case of a spring, 


8(x - x-l) a function such that its value is zero at every point except 

the argument x^ and its value at this point is infinity in 

I* s(x - Xj_)dx =1 

-L 


such a way that 


The coupled bending and torsional deformations y and cp can be 
expressed in terms of the coupled modes of the unmodified beam by 
infinite series 


7 = WO + Wl + • • 


<P = + ^^l + ' * 

■where the coefficients a^, a 1 , . . . 

coordinates that give the proportions of the primary-beam modes that are 
present in a natural mode of the modified structure. 

Since y m and cp m are coupled natural modes of the primary 
structure, ,they must satisfy the differential equations 



and the orthogonality conditions, if m £ n. 


00 

= 2 


ut=0 

OO 


■I 


nt=0 


Smyin 


a cp 
nr m 


( 3 ) 


can be considered as generalized 


J + e<pn)dx + J mcp m (k?cp n + eyjdx =0 (5a) 

— I* — Xi 
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and, if m = n. 



In addition, the pr i ma r y- structure modes y n and cp satisfy the same 
boundary conditions as the modified-structure modes y and cp. 

The procedure -which is used to obtain the characteristic equation 
is based on the Gale r kin method and is as follows: The infinite series 

for y and cp (equations (3)) are substituted into the differential 
equations (la) and (lb); then equations (la) and (lb) are multiplied by 
y n and cp n , respectively, and are integrated over the length of the 
beam. This procedure yields 



Ii 00 

J yn™^ ^frm + ec Pm) dx 

“-k m=0 




y n C p ]i 6 C x - x ±)d* = 


a 


(n = 0, 1, 2, . . . ) (6a) 


and 




^ ^(^m + k^dx 


m=0 


/ f 'Pn&Ij. 6 ^ " = 0 (n = 0, 1, 2, . . .) (6b) 

i ^ 


With the use of the definition of 5(x - x-j_) , the differential 

equations shown in equations (4), and the orthogonality conditions in 
equations ( 5 ), the subtraction of equation (6b) from (6a) can be 
simplified to 


an^nta 2 - <0 - ± [l n ] ± + Mi^n]^) (7) 

(n = 0, 1, 2, . . . ) 
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If [p^ and Mi are expanded in terms of a n y n and a^cp^ 

(equations (2) and (3)), equations (7) would lead to a set of simultaneous 
homogeneous equations in the coefficients a^. The order of the frequency 

determinant that would result from the direct consideration of this set 
of equations would he equal to the number of terms used in the deflection 
functions (equations (3))- A simplified frequency determinant can be 
developed, however, in the following manner. 

Simplified frequency determinant .- The solution of equation (7) 
for a^ gives 


station x^ can be written 


The loads and w due to any modification at a particular 

CO 

Eh = XI ^Pnl, (j = 1, 2, . . . 1) (9a) 

J n=0 L 


and 


CO 

t], = 2Z anfn], (J = 1, 2, . . . i) ( 9 b) 


n=0 


where 


" K i (M 


j +e J 


Nj = + k J 


The substitution of equation (8) into equations ( 9 ) gives 


Mj) 
Mj) 


(10) 


Mj - 51 (?], 2 




n=0 rf n (o^ 2 - i 


SI 5 


a=0 N^ 2 - a £) 

[ Tfl ]jWi 


1 i n=0 N n (c% 2 - of) i n=0 N^o^ 2 - cn 2 ) 


(U) 


where (p] ^ and pfj -j_ are unknowns . 


(J = 1, 2, . . . i)' 
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Equations (ll) can be written for every station at which a modifica- 
tion is made; thei-efore two equations are obtained for every concentrated 
mass or spring on the structure. These equations are shown in table 1. 

Since these equations are homogeneous in. [p^ and Qr]^, the 

determinant of these equations, which is the bracketed part of table 1, 
must vanish in order to have real values of DGi and DOf* The order 

of this simplified frequency determinant is now proportional to the 
number of modifications instead of the number of terms in the expansion. 
Although each term of the determinant is an infinite series, for practical 
problems it is necessary to include only the first few terms of each 
series, each term of which corresponds to a known mode of the unmodified 
structure. As the examples show subsequently, at least as many frequencies 
and modes of the unmodified beam must be known as the number of frequencies 
and modes of the modified beam desired. The rapid convergence of the 
calculated frequencies to the exact frequencies as the number of terms of 
the series is increased is shown in the examples. 

Calculation of the modes and frequencies .- For most cases perhaps 
the simplest method of obtaining the natural frequency cn of the modi- 
fied beam from the frequency determinant is to evaluate the determinant 
for several trial values of ■ a) in the expected vicinity of the natural 
frequency and plot the value of the determinant against the trial fre- 
quency cs. In most cases only three or four points are needed to obtain 
the natural frequency. 

This method of evaluating the natural frequencies is especially 
advantageous when the effect of several loading conditions in which the 
modifications do not change position but only magnitude is being studied. 
Under this condition the value of each summation in all the elements of 
the determinant change only in proportion to the changes in the values 
of K ± . 

To calculate the natural mode after its frequency has been established, 
the relative values of and (tQ^ have to be obtained. In order to 

find these relative values, one of the loads is given an arbitrary value 
of unity and the resulting set of north omogeneous equations is solved 
simultaneously for the relative values of the other loads. In the process 
one of the equations can be discarded, but it has been found best simply 
to add two of the equations together before the solution is made. 

With the values of the frequency and the loads [p] ± and (50^, 

equation (8) is used to evaluate the coefficients Sq, and then the modal 

functions can be found from equations (3)* At least the same number of 
primary-beam modes should be included in the calculation of the modified- 
beam mode as terms used in the evaluation of the series. 
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The convergence of the calculated mode to the exact mode as the 
number of primary-beam, modes included in the calculation is increased 
is shown in the examples. 


EXAMPLES 


In order to illustrate the method and to show the accuracy with 
which the frequencies and modes of the modified beams can be calculated, 
the uncoupled symmetrical bending frequencies and modes of a uniform 
free-free beam are calculated with the following modifications made at 
the elastic axis of the beam (see fig. 3): 

Example 1 - Equal masses placed at the tips of the beam 

Example 2 - Masses placed at the center line of the span and third 
points of the beam 

Example 3 - A spring placed at the center line of the beam 

Because the masses and springs were placed at the elastic axis, e i 
and are zero or 

Mi = K i[y]i (12a) 

and 

Hi = 0 (12b) 


Also, because only the symmetrical uncoupled modes are required, only 
symmetrical uncoupled modes and frequencies of the primary beam are used. 

Modes and frequencies of the unmodified beam .- The frequencies of 
the unm odified uniform free- free beam are 



, k EI 

11 mT. ^ 


f^ = 31.28 

f 2 k = 913 . 6 


f 3 k = 5,571 

f k k = 19,262 

f^ = b- 9,600 


where 



10 


NACA TN 2132 


The constants N are 


N Q = 2mL 


and, for n ^ 0, 



The modes of the unmodified uniform "beam, normalized to give a unit 
tip deflection, are given by the following equations: 


^n o l 


/cosh f n - cos f n f 
L + L 


cosh f„ 


cos f. 


n 


(n = 0, 1, 2, . . . ) 


Example 1: Calculation of the frequencies and modes of a uniform 

free-ft~ee beam with concentrated masses at the tips (fig. 3(h)) .- 
Example 1 illustrates the use of the method- for the case of symmetrical 
modifications. The placement of the masses at the tips gives a severe 
test of the accuracy of the method. The results obtained for this example 
are compared with an exact solution of the differential equation. 

Since the modifications are symmetrically placed, that is, x-^ ■= - L , 

Mi = M a 

With the use of this condition and equation (12), the frequency determinant 
for equal tip masses can be written as 



The results obtained from the solution of this frequency equation when 

tnT 

— = ^ are given in the following table for 


CD 


mL 
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Mode 

Number of terms 
in series 

Calculated f^ 

Exact f 2 

Percent 

error 

First 

3 

4.15 

4.15 

0 

First 

2 

4.17 

4.15 

• 5 

Second 


25.0 

24.9 

.4 

Second 

3 

25-2 

24.9 

1.2 

Third 

6 

65.0 

64.8 

.3 

Third 

5 

65.2 

64.8 

.6 

Third 

. * 4 

65.8 

64.8 

1.5 


Examination of this table shows that the use of the series with a 
mi nlmum number of terms evaluates the frequency of first three modes 
with errors of less than 2 percent, which is sufficient for most uses. 

Once the frequencies of the modified beam are found, the proportion 
'of each of the primary-beam modes present in a modified-beam mode is 
calculated from equation (8). Because the values of [p]^ are equal, 

equation (8) becomes 


N n (<% 2 - <*&) 


after M i is assumed to be unity. 

The comparison of the solution for the first three modified-beam 
modes with the exact modes is shown in figure 4. 

Example 2: Calculation of the frequencies and modes of a uniform 

free-free beam with concentrated masses at the center line and third 
points of the beam (fig. 3(c) ) «- Example 2 illustrates the use of the 
method when more than one modification is made to the primary structure. 
The results for this example are compared with calculations obtained 
from a matrix iteration process which is considered as exact. 

The magnitudes of the masses and the assumed geometrical properties 
of the beam are shown in figure 3(c)* Since, as in example 1, two of 
the masses are symmetrically placed, the frequency determinant can be 
reduced to 
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Mi 


-1 + 2 




n=0 


N 


*c[?n ] C 2 

n(<% 2 - <^) 


5 




n=0 »n(% 2 - m 2 ) 


Ha 


-y- ^ggAgjc 

^CT ®n(<% 2 - °?) 


= 0 


-1 


+ 


00 


2 


%| yJ A 2 

u n («% 2 - “ 2 ) 


The results obtained from the solutions of this frequency determinant 
are as follows: 


Mode 

Number of terms 
in series 

Calculated o> 

'Exact co 

Percent 

error 

First 

3 

21.8 

21.8 

0 

First 

2 

22.0 

21.8 

•9 

Second 

5 

91.4 

91.4 

0 

Second 

4 

93-0 

91.4 

1.7 

Second 

3 

100.7 

91.4 

10.2 

Third 

5 

261.8 

261.3 

.2 

Third 

4 

269.5 

261.3 

3-1 


For this example equation (8) becomes 

«■> --rrt — + mawO 

N n (t% 2 " ® < V N 

when ^ is ass ume d to be- unity and the relative value of Ha 18 

calculated from the deter minan t just shown, in the manner previously 
described. The comparison of the calculated and the exact solutions for 
the first three modes is shown in figure 5. 

The large error in the calculation of the frequency of the second 
mode with the minimum number of terms is caused by the large peak 
deflection that occurs in the second mode near. the point where the masses 
are placed. 
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Example 3: Calculation of the frequencies and modes of a uniform 
free-free beam with a spring at the center line .- Example 3 is presented 
to show how the method is used when a beam is modified by the addition 
of a spring instead of a mass. The results obtained are compared with 
an exact solution of the differential equation. 

Because only one modification is used, the frequency determinant 
consists of only one term 


1 + 



0 


The results obtained from the solution of this frequency equation when 
a spring is used with a spring constant 



are tabulated as follows for 


oi 2 = f^ ~ 

mL 


4 


Mode 

Number of terms 
in series 

Calculated f2 

Exact f2 

Percent 

error 

First 

3 

2.49 

2.49 

0 

First 

2 

2.51 

2.49 

.8 

Second 

4 

7.01 

7.00 

.1 

Second 

3 

7.02 

7.00 

• 3 

Third 

5 

30'. 6 

30.6 

0 

Third 

4 

30.6 

30.6 

0 


The proportion of each of the primary-beam modes present in the 
modified-beam modes, as in other examples, is obtained from equation (8), 
which can be written 


B n=~ 
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The comparison of the solutions for the first three modi fled- he am 
modes with the exact solutions is shown in figure 6. 


CONCLUDING REMARKS 


The method presented herein can he used once the frequencies and 
modes of the unmodified structure are known. Thus the method is 
advantageous when the modes and frequencies of the unmodified structure 
either are known or are required along with the modes and frequencies 
of the structure with various modifications. 

The amount of time necessary to calculate a natural frequency and 
mode of a modified structure once the primary- structure modes and 
frequencies are known by the method presented herein is generally shorter 
than any of the standard methods. Also, in most cases that can he 
encountered, all the calculations necessary can he done with a slide rule. 
Most nf the calculations done for one frequency and mode can he used in 
the calculations of other frequencies and modes. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , April 2k, I 95 O 
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Figure I. — Structure used in theoretical analysis. 



Figure 2. — Definition of the eccentricities. 
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* 2 

L > 


_^x 


<L 

(a) Unmodified uniform free-free beam. 
M M=0.25mL M 



(b) Example I — Addition of masses at tip. 



02 5 M A = M b = Mcf 0.9 42 

TlS ' 5 - 514 

(c) Example 2-Addition of multiple masses. 



777777T 

(d) Example 3 — Addition of a spring. 


Figure 3.— Unmodified beam with modifications used 
in examples. 



Relative deflection Relative deflection 
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Figure 4.— Comparison between calculated and exact solutions in example I. 
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Fiqure 5. -Concluded 





Relative deflection Relative deflection 
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Exact solution 
Solution using 2 modes 



(a) First mode. 


Exact solution 
Solution using 3 modes 


(b) Second mode. 


Figure 6. — Comparison between calculated and exact solutions in example 3 
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Figure 6.- Concluded. 
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